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ABSTRACT 

The first detection of the Cqq (Buckminsterfullerene) molecule in massive embed- 
ded young stellar objects (YSOs) is reported. Observations with Spitzer IRS re- 
veal the presence of C 60 in YSOs ISOGAL-P J174639.6-284126 and SSTGC 372630 
in the Central Molecular Zone in the Galactic centre, and in a YSO candidate, 
2MASS J06314796+0419381, in the Rosette nebula. The first detection of C 60 in a 
Herbig Ac/Be star, HD 97300, is also reported. These observations extend the range 
of astrophysical environments in which Cgo is found to YSOs and a pre-main sequence 
star. Ceo excitation and formation mechanisms are discussed in the context of these 
results, together with its presence and processes in post-AGB objects such as HR 4049. 

Key words: astrochemistry - circumstellar matter - Galaxy: centre - ISM: molecules 
- stars: pre-main-sequence - stars: winds, outflows 



1 INTRODUCTION 

Although considered a candidate carrier for a few infrared 
astrono mical emission f eatures, as noted in Werner et alj 
(2004b) and discussed bv lSellgren. Uchida fc Wernerl l|2007f) . 
confirmation o f Ceo together w ith C70 in the planetary 
nebula Tc 1 |Cami et al.l bold ) and C60 in the reflec- 
tion n ebulae NGC 7023 and NGC 2023 (|Sellgren et all 
l20ld . l201ll ) has only recently been achieved. Fur- 
ther discoveries of neutral Ceo have since been re- 
ported in planetary nebul ae (PNe) in the Milky Way 
and the Magellanic Clouds l|Garcfa- Hernandez et al.ll2010l. 
l201ll). weakly H-deficient R Coronae Borealis- t ype s tars 
( Garcia- Hernandez. Kameswara Rao &: Lambertl |201ll ) . a 
protoplanetary nebula ( Zhang fc Kwokl 2011), ac ross the 
'veil' region of the Orio n Nebula |Rubin et a l *l 20 111 ), in sev- 
eral post-AGB objects dGielen et al.l|2011al lbh and possibly 
in the binary XX Oph (jEvans et al.ll201ll ). There is there- 
fore now unambigious evidence for neutral C60 in objects in 
the later stages of stellar evolution and in the interstellar 
medium. 

We report here the first detection of C60 in two young 
stellar objects (YSOs) located in the Central Molecular Zone 
(CMZ) of the Milky Way, in a YSO candidate in the Rosette 
nebula and in a Herbig Ae/Be star. Taken with earlier re- 
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ports, these findings show that neutral C60 exists over a very 
wide spectrum of stellar and interstellar evolution, encom- 
passing star forming regions and young stars, mass-losing 
evolved stars, protoplanetary and planetary nebulae, and 
the interstellar medium. We also report detection of C60 in 
the post-AGB star HR 4049 which , like HD 52961 in which 
C60 is known l|Gielen et al.ll2011bl ). is unusual in displaying 
the rare infrared emission signatures of nanodiamonds. 

This paper is organised as follows: background to cur- 
rent issues associated with astrophysical C60 is presented in 
section [2] details of the Spitzer observations and data re- 
duction in section[3]and the classification of the YSO£| and 
the YSO candidate in section f5] Measurements and analysis 
of C60 features are presented in section [5] Excitation mech- 
anisms and C60 formation are discussed in section [6] and the 
results are summarised in section [7] 



2 BACKGROUND 

Astronomical conditions analogous to those of t he early lab- 
oratory experiments in whic h C60 is formed l|Kroto et alj 
ll985l ; lKratschmer et al.lll990h might have been expected to 
favour C60 formation, with hydrogen-poor R Coronae Bore- 
alis (RCB) stars being good candi dates. However, this has 
not been borne out by observation (|Clavton et al.lll995l ) ex- 
cept in the case of one or possibly two less H-poor RCB stars 



1 We choose to use the term YSO to describe the embedded YSOs 
rather than include the more evolved Herbig Ae/Be star. 
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(IGarc i'a-Hcrna ndez. Kameswara Rao fc Lambert! lioill ). Al- 
ternatively, laser-induced decomposition of hydrogenated 
amorphous carbon (HAC) i s known to produce fullerenes 
|Scott. Dulev fc Pinholll997l ') and HAC decomposition has 
been r eferred to as a poss i ble astrophysical Ceo formation 
route (|Moutou et al. 1 1 19991; lEhrenfreund fc Charnlevl I2OO0I; 



ISellgren et al. 20101 ; Garcfa-Hernandez et al. 2010T ). Follow- 
ing identification of C60 and C70 in the particularly clean 
spectrum of Tc 1, ICami et al.l (|2010h suggested that these 
molecules formed efficiently in this object because the cir- 
cumstellar environment was H-poor. However, in their pa - 
per on C60 emission in PNe[G arcia- Hernandez et all (<2010l l 
pointed out that neither the PN nor its compact core and 
the central star of Tc 1 are H-poor, and presented evi- 
dence in support of the idea that fullerenes are formed 
in astrophysical environents by decomposition of HAC 
|Garcfa-Hernandez et al.l l201lT ). They suggested that C60 
and polycyclic aromatic hydrocarbons (PAHs) were likely 
to be formed together from the decomposition of HAC due 
to UV processing and energetic phenomena such as shocks, 
and interpreted the deficiency of PAHs in Tc 1 in terms of 
longer survival times of fullerenes due to their high stabil- 
ity. C60 formation issu es are further explored in papers by 
ICami et al.l l|201ll ) and ICarcfa-HernandezI (|201lh . with the 
challenges being extended further by the detection of C 60 in 
mixed-chemistry post-AGB stars (I Gielen et al.ll2011al ) and 
reported for HR 4049 in this paper. 

The phase (gas or solid) and excitation mechanism of 
Ceo in astrophysical sources is not firmly established and 
may vary between environments. From their analysis of spec- 
tral data for Tc 1 and particularly the relativel y low vibra- 
tional temperature of ~ 330 K that was inferred. ICami et alj 
(2010) suggested that C60 in Tc 1 is attached to, and in 
thermal equilibrium with , solid materia l such as carbona- 
ceous grains. In contrast. ISellgren et al.l (|2010T l have found 
the vibrational band intensity ratios in NGC 7023 to be 
consistent with high-energy (~ 10 eV) UV photo- excitation 
of gas-phase Ceo- At present we are not aware of predicted 
band strengths based on other possible excitation mecha- 
nisms. The discovery of C60 in YSOs and a Herbig Ae/Be 
star provides additional environments with which proposed 
formation and excitation mechanisms can be evaluated. 



3 OBSERVATIONS AND DATA REDUCTION 

A search was undertaken for C60 emission bands 
in sp ectra recorded w ith the Infra-Red Spectrograph 



(IRS; iHouck et al.l |2004 ) of the Spitzer Space Telescope 
l| Werner et al.ll2004al ). We interrogated the Spitzer Heritage 
Archive (SHAjf for objects observed at the wavelengths of 
the C60 emission bands, concentrating on pre-main sequence 
objects and examining several hundred spectra. Targets were 
selected by visual inspection of pipeline-reduced data in the 
spectral region around the 18.9 C60 band, which is rel- 
atively strong and free from blending with atomic lines and 
PAH emission bands. Data for each object with a potential 
18.9 (ixn detection were retrieved from the archive, together 
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with additional targets selected for comparison. Observa- 
tional details and photometry for the objects found to con- 
tain C60 are given in table [1] 

All of the targets were observed using the IRS Short- 
High (SH) module which covers the wavelength region 9.9- 
19.6 jrni with resolving power R = ~ 600. Some targets 
were also observed with the Short-Low (SL) module (5.2- 
14.5 jmi, R ~ 60-120), and/or the Long-High (LH) module 
(18.7-37.2 jim, #~ 600). 

Flux- and wavelength-calibrated Post Basic Calibrated 
Data produced by version S18.18.0 of the Spitzer Science 
Centre pipeline were obtained from the archive, with fur- 
ther analysis undertaken using iRAljf]. The observations were 
backgro und-subtracted after inspection of 2MASS K-hanA 
images jSkrutskie et alj|2006j ) for the presence of contami- 
nating sources. Where no such sources were present, either 
the background subtraction performed by the pipeline was 
used (SL), or all of the background regions were combined 
and manual subtraction was employed (SH and LH). Where 
such sources were found to be present, only those back- 
ground regions without significant contamination were com- 
bined and the background was subtracted manually. Multi- 
ple exposures were combined to increase the signal-to-noise 
ratio. 

For targets observed with more than one IRS module, 
the spectra were combined by scaling the flux in overlap- 
ping wavelength regions. Where these were in disagreement, 
the flux from the SH observations was preferred and the 
other fluxes scaled accordingly, because the SH module has 
a smaller field of view and thus samples the flux from the 
target more reliably. In the final spectra, the wavelength 
ranges used from each module were: 5.2-10 (SL), 10-19.5 
(SH) and 19.5-36 jrni (LH). 



4 OBJECT CLASSIFICATION 

Given the importance of determining that the YSOs dis- 
cussed in this paper are embedded young objects, we re- 
view here the evidence that our two CMZ objects, ISOGAL- 
P J174639.6-284126 and SSTGC 372630, are indeed YSOs 
and not dust-enshrouded post-AGB sources. Based on the 
available data, we designate our third (Rosette nebula) ob- 
ject, 2MASS J06314796+0419381, as a candidate (but very 
likely) YSO. The attribution to YSOs rests principally on 
the steep smooth rise in the spectral energy distribution 
(SED) with increasing wavelength, the presence of silicate 
and CO2 ice absorption bands, the [Nell] 12.8 Jim line and 
PAHs in emission, and the absence of a broad 30 (im emis- 
sion featurjfl which is well correlated with the presence of 
C 6 o in post-AGB objects. The Herbig Ae/Be star HD 97300 



3 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

4 The 30 urn f e ature has been ascribed to MgS 
llGoebel fc Moselevl 1 19851 ; iHonv, Waters fc Tielensl 120021) but 
may, perhaps more probably, be due to carbon aceous material 
llZhang, Jiang fc Lil I2OO9I ; IZhang fc Kwokl I201ll 1 such as HAC 
iGrishko et alj|200ll) . 
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Table 1. Coordinates and photometr y for the targets with C 60 emission bands. Right ascension, declination, and near-infrared magnitudes 
are taken from the 2M ASS catalogue ( Skrutskic ct aL 2006); mid-infrared magnitudes (where available) are from the Spitzer IRAC survey 
of the Galactic centre toarm'rez et alT T2008). Spitzer programme numbers and Principal Investigator (PI) names are are given for the 
IRS observations used in this study. 



Name 



RA 



Dec 



A' 



[3.6] [4.5] [5.8] [8.0] Programme (PI) 



Embedded YSOs: 
ISOGAL-P J174639.6-284126 
SSTGC 372630 
2MASS J06314796+0419381 



17:46:39.60 -28:41:27.0 > 13.8 12.95 10.28 8.83 7.38 5.58 
17:44:42.76 -29:23:16.2 > 16.0 12.87 10.31 8.82 7.67 6.48 
06:31:47.96 +04:19:38.2 14.01 10.67 - 



40230 (Ramirez) 
40230 (Ramirez) 
50146 (Keane) 



Other targets: 

HD 97300 (Herbig Ae/Be star) 11:09:50.03 -76:36:47.7 7.64 7.15 

HD 52961 (post-AGB object) 10:18:07.52 -28:29:30.7 16.06 15.42 

HR 4049 (post-AGB object) 07:03:39.63 +10:46:13.1 6.32 5.53 



2 (Houck) 

3274 (Van Winckel) 
93 (Cruikshank) 



and the post-AGB star HR 4049 are well known and not 
discussed in this section. 



4.1 ISOGAL-P J174639. 6-284126 

ISOGAL-P J174639. 6-284126 is located in the CMZ of the 
M ilky Way (see tab le [P. It was identified as a massive YSO 
bv lFelli et all (2002) based on ISO mid-infrared photometry. 
Its Spitzer IRS spectrum is shown in figure[T]as trace (a) and 
illustrates a steep rise in SED with increasing wavelength, 
an SED peaking to the red of ~35 um being indicative of 
a young object (e.g. ISimpson et al.ll201ll ). The figure shows 
that the SED of ISOGAL-P J174639.6-284126 has a much 
faster rise at longer wavelengths than those of t wo sample 
post-A GB objects shown in traces (c) and (d). lAn et al.l 
i|201ll) identified ISOGAL-P J174639. 6-284126 as a candi- 
date YSO based on 2MASS and Spitzer photometry, but 
excluded it from their final list of YSOs due to the lack of 
a ~ 15.4 urn shoulder on the CO2 ice absorption feature. 
The presence of a shoulder at ~ 15.4 urn, thought to be 
due to methanol-C02 ice complexes, was used as the YSO 
selection criterion principally to distinguish between YSOs 
and field s tars behind molecular clouds. Howeve r, other au- 
thors (e.g. lSeale et a l. 2009; Simp son et al.ll201ll ) have iden- 
tified massive YSOs which do not have CO2 ice absorp- 
tion features, or have different peak w avelengths and pro- 
files f rom those used by An et al.l (2011) for selecting YSOs. 
Also, lEhrenfreund et aL ( 19991 ) have shown that photoly- 
sis of methanol-containing CO2 ice mixtures decreases the 
strength of the shoulder on the CO2 ice band; this could lead 
to disappearance of the shoulder from spectra of the more 
evolved YSOs, such as those associated with a compact or 
ultra-compact Hn region that are likely to have higher UV 
flux (see section l4~4)l . 

From figure [1] it is clear that there is no 30 um 
emission feature a s seen in the Cgo-containing PPN 
IRAS 01005 +7910 (IZhang fc Kwokl [20111) and PN SMP 
SMC 016 (|Garcfa-Hernandez et all l20ld )Fl Figure [2] 
shows the continuum-normalised spectrum of ISOGAL- 
P J174639. 6-284126 between 10 and 19.5 um with two sil- 
icate absorption bands and absorption by CO2 ice indi- 



SEDs of YSOs and dusiy posi-AGB objects 
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5 Spectra of P Ne llBernard-Salas et alj |2009|) and PPNe 
dVolk et al.ll201ll) provide further examples. 



Figure 1. Spitzer IRS spectra of the CMZ objects (a) ISOGAL-P 
J174639. 6-284126 and (b) SSTGC 372630, illustrating the strong 
rise in the far-IR SED which is indicative of a YSO. The compar- 
ison traces are for (c) proto-planetary nebula IRAS 01005+7910 
and (d) planetary nebula SMP SMC 016. These show a slower rise 
in SED (indicated by the dashed lines) and a broad 30 (J.m emis- 
sion feature commonly seen in C60-containing post-AGB objects. 
The spectra are normalised to the flux of ISOGAL-P J174639.6- 
284126 at 22.5 [im and vertical offsets applied. 



cated. Based on these spectral characteristics, together with 
the presence of PAH and [Nell] emission in this (and our 
other) YSOs which is generally absen t in low-mass YSOs 
|van Dishoeckll2004l ; iGeers et alj|2009l ). it is concluded that 
ISOGAL-P J174639.6-284126 is a high-mass YSO. 



4.2 SSTGC 372630 

SSTGC 372630 is also located in the CMZ (see table Q). 
Based on its SED and position in an infrared colour- 
magnitude diagram, it w as identified as a massive YSO by 
lYusef-Zadeh et al.l (|2009r ). Its Spitzer IRS spectrum between 
5 and 36 urn is given in figure[T]- trace (b) and, like ISOGAL- 
P J174 639. 6-284126, sh ows a steeply rising continuum in the 
far-IR. lAn et all |201ll ) listed this object as a possible YSO 
based on photometric data and their analysis of the shape 
of its strong CO2 ice feature at ~ 15.2 um. Our own data re- 
duction does suggest an absorption component near 15.4 um. 
There is very strong absorption due to silicates, CO2 ice at 
~15 um, ices at 6.0 and 6.8 um and PAH emission. There is 
also H2 rotational line emission which is an indicator of an 
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Figure 2. Continuum-normalised spectrum of ISOGAL-P 
J174639. 6-284126 between 10 and 19.5 um showing the presence 
of silicate and CO2 ice absorption. 



2MASS 106314796+0419381 




Figure 3. Continuum-normalised spectrum of 2MASS 
J06314796+0419381 between 10 and 19.5 um showing the 
presence of silicate and CO2 ice absorption. 



embedded YSO which has not yet dev eloped an ultracom- 
pact Hll region (|Varricatt et al.ll201Cf ). It appears to have 
th e spectral chara c terist ics of an 'outflow' YSO as described 
bv lSimpson et all (|201ll ). We conclude that SSTGC 372630 
is also a massive YSO. 



4.3 2MASS J06314796+0419381 

2MASS J06314796+0419381 lies in the outskirts of the 
Rosette nebula, which is a site of recent star formation 
IlKuchar fc Baniall 19931 ). Invisible in optical images, this ob- 
ject only becomes detectable in the near-IR (see tableQJ and 
brightens rapidly longward of 15 urn. The heavy extinction 
of this object, together with CO2 ice and silicate absorption 
(see figure [3}, the presence of H2 rotational emission and 
PAH features are supportive of it being classified as a mas- 
sive YSO. There are no photometric or spectroscopic data 
available beyond 21 urn so the presence or absence of a 30 um 
emission feature is not known. As there are no directly rel- 
evant photometric data and the longer wavelength SED is 
not established, we treat 2MASS J06314796+0419381 as a 
candidate massive YSO. 



4.4 YSO Evolutionary stage 

The three YSO targets are at different stages of evo- 
lution, and fall between the extremes represented by 
Cep A East and S106 IR which are examples of deeply 
embedded and more evolved massive YS Os, respectively 
l|van den Ancker. Tielens fe Wesselius|[2000T ). In Cep A East 
there are strong silicate and ice absorption features but no 
PAHs, whereas S106 IR is characterised by emission in fine- 
structure lines of atomic ions and from PAHs, reflecting the 
impact of UV radiation on the environment and its capacity 
to create Hn regions and PDRs. Emission from atoms and 
molecules excited through shocks of various types is also 
expected to be prevalent. In broad terms SSTGC 372630 
is relatively deeply embedded, 2MASS J06314796+0419381 
is less so and ISOGAL-P J174639.6-284126 lies closer in 
form to the more evolved object S106 IR. The more evolved 
stage of ISOGAL-P J174639.6-284126 is supported through 
prominent [S ill] (very weak in our other two objects), strong 
[Ne 11] and PAH band ratios similar to those seen in H n re- 
gions. It is likely that both shocks and UV flux are present 
in the three young objects but to differing degrees. 



5 SPECTRA OF C 6 ,,-CONTAINING OBJECTS 

In this section we describe analysis of the C60 emission bands 
in the two YSOs and candidate YSO, the Herbig Ae/Be star 
HD 97300 and the post-AGB object HR 4049, spectra of 
which are shown in figures 4-8. For each emission feature, 
integrated fluxes (table [2} were measured after subtraction 
of the nearby continuum which was fitted with a cubic spline. 
Whilst the C60 band at 18.9 um is an isolated feature, the 
other three C60 bands are blended with atomic and/or PAH 
emission to varying degrees for each object. The 7.0 um Ceo 
band is sometimes contaminated by an [Aril] emission line 
at 6.99 um and potentially also by a band of C70, but as 
this region is only available in the low resolution SL spectra 
deblending of these features is difficult. The 8.5 um C60 band 
is heavily blended with a broad PAH feature peaking at 
8.6 urn which has prevented measurement of the flux in this 
band for all targets. The 17.4 um C60 band may include some 
contribution from PAH emission. None of the spectra have 
the strongest C70 emission feature at 15.6 um, suggesting 
that C70 is either absent or present in concentrations too 
low to be d etected, and no evidenc e for bands of C24 as 
predicted bv lKuzmin fc Dulevl (|201ll ) was found. 

The YSOs with Ceo bands (sections l5.H5.3p have a 
number of common spectral characteristics. All have PAH 
features and [Nell] line emission, and two of the objects 
(SSTGC 372630 and 2MASS J06314796+0419381) have 
pure rotational H2 emission lines. In ISOGAL-P J174639.6- 
284126 there are also [Sin] (34.82 um) and [S ill] (18.71 & 
33.48 um) emission lines. 



5.1 YSO ISOGAL-P J174639. 6-284126 

The spectrum of ISOGAL-P J174639. 6-284126 covering the 
Ceo bands is shown in figure [4] and the derived fluxes are 
listed in table [5] In this and for all other objects the flux 
quoted includes any [Aril] contribution to the 7.0 um band. 
The absence of the high ionisation lines - [Neni] 15.55, 
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Table 2. Wavelength (A, in (im), full width at half maximum (FWHM, in urn) and flux (in mjy urn) for the 7.0, 17.4 and 18.9 um C60 
emission bands in the YSO and other targets. 



Object 




7.0 um 




17.4 um 




18.9 urn 




A 


FWHM 


Flux / mjy um 


A 


FWHM 


Flux / mjy um 


A 


FWHM 


Flux / mjy um 


ISOGAL-P J174639. 6-284126 


7.01" 


0.09 a 


57.7 ± 1.0 a 


17.38 


0.29 


72.1 ±8.5 


18.91 


0.32 


169.5 ±7.9 


SSTGC 372630 


7.04 


0.16 


13.28 ± 1.1 


17.38 


0.34 


12.9 ±1.3 


18.92 


0.30 


21.6 ± 1.2 


2MASS J06314796+0419381 


7.01 


0.11 


6.47 ±0.14 


17.38 


0.22 


40.550 ± 1.9 b 


18.94 


0.30 


61.7 ± 2.7 


HD 97300 














18.92 


0.27 


61.8 ± 1.7 


HR 4049 








17.32 


<1.24 


< 1239 ± 28 


18.89 


0.64 


1565 ± 48 


HD 52961 c 








17.39 


0.28 


37.5 ±2.2 


18.91 


0.35 


99.0 ± 1.9 



Includes a contribution from an [Aril] line. b Includes a contribution from a PAH feature. c See also Gielene^aL 



ISOGAL-P .11 74f-o l i.(>-iS4i:(i 



SSTGC 372630 




Wavelength / Jim 

Figure 4. Spitzer IRS spectrum of the CMZ YSO ISOGAL-P 
J174639.6-284126 cov ering the region o f the C 6 o bands at 7.0, 
8.5, 17.4 and 18.9 um llFrum et al.|jl99lt) . The expected positions 
of the C60 bands are marked with arrows, which are solid where 
the feature is detected and dashed where only an upper limit on 
the flux is determined. PAH features at 6.2, 7.7, 8.6, 11.2, 12.7 
and 16.4 um are also marked. 
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Figure 5. Spitzer IRS spectrum of the CMZ YSO SSTGC 372630 
covering the region of the C60 bands. In addition to the lines 
identified in figure[4] very strong absorption due to silicates, CO2 
ice at ~15 um and ices at 6.0 and 6.8 um is present. Rotational 
lines of H2 are also marked. Note the very strong PAH emission 
band at 7.7 urn. 



[Arm] 8.99, and [P hi] 17.89 u.m - whic h are seen in the 
planetary nebula Tc 1 (jCami et al.ll2010h suggests that the 
UV radiation field is less harsh than in Tc 1. However, 
the observed [Nell] 12.81 and [S ill] 18.71 urn lines do indi- 
cate the presence of a compact or ultracompact H n region 
ijSeale et al.ll2009l ). The PAH spectrum, and in particular 
the low 11.2 (j.m/12.7 u.m band ratio is also con sis- 

tent with the presence of a H 11 region (|Honv et al.ll200lh . 

5.2 YSO SSTGC 372630 

In addition to Ceo emission in the 18.9 and 17.4 (im bands, 
the spectrum of SSTGC 372630 contains Ceo emission at 
7.0 urn, strong silicate, CO2 and other ice absorption fea- 
tures, and H2 S(l) and S(2) rotational emission lines (see 
figure[5]). The 7.7 urn PAH emission band is distorted by the 
broad silicate absorption, but still appears very strong when 
compared to the other PAH bands. Although the [Nell] line 
is much weaker than in ISOGAL-P J174639.6-284126, there 
may still be a small contribution to the 7.0 urn emission fea- 
ture from [Aril]. 

5.3 YSO Candidate 2MASS J06314796+0419381 

The spectrum of 2MASS J06314796+0419381 given in fig- 
ure [6] has the weakest ionised atomic lines of the three ob- 



2MASSJ063 14796+0419381 




6 8 10 12 14 16 18 

Wavelength / fun 

Figure 6. Spitzer IRS spectrum of 2MASS J06314796±0419381 
in the Rosette nebula covering the region of the Cgo bands. An- 
notations are as in earlier figures. There is a rapid rise in flux at 
longer wavelengths. 



jects, with barely discernible [S ill] , weak [Ne n] and hence 
likely negligible [Am] at 6.99 urn. Strong Ceo bands at 18.9, 
17.4 and 7.0 u.m are readily identifiable, together with PAH 
emission bands and the S(l), S(2) and S(3) lines of H2, where 
the S(3) line falls in the region of very strong silicate absorp- 
tion. 
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HR4049 and HD ?2 ( )hl 




Figure 7. Spitzer IRS spectrum of the Herbig Ae/Be star 
HD 97300 covering the region of the C60 bands. Annotations are 
as in earler figures. 




Figure 8. Spitzer IRS spectra of the post-AGB objects HR 4049 
and HD 52961. The * symbol denotes gas-phase CO2 emission. 
Inset: a X 10, offset 60 and b X 100, offset 110. 



5.4 Herbig Ae/Be star HD 97300 

HD 97300 is an optically bright Herbig Ae/Be star of 
spectral type B9. It is a more evolved system than 
the YSOs and is approa ching the stellar main sequence 
|Siebenmorgen et al.lll99a ). The IRS spectrum is shown in 
figure [3 where C60 emission in the 18.9 um band is visible. 
Emission is also seen in the region of the 17.4 u.m Ceo band 
but this is heavily contaminated by PAH bands, so it is not 
possible to measure the flux from Ceo- No C60 feature is 
seen at 7.0 u.m suggesting a low level of vibrational excita- 
tion, correspondin g to a vibrational temperature of 200 K 
l|Cami et al.l 1201 ll ) - see also section 16.11 Weak emission in 
the S(l), S(2) and S(3) lines of H2 is present. Although the 
identification of Ceo in HD 97300 is based mostly on the 
18.9 u.m emission feature, it is reproducible as seen in other 
spectra dKeller et alj|2008l ; iManoi et alj|201ll ). 



5.5 Post-AGB star HR 4049 

The Spitzer spectrum for this post-AGB object is shown in 
figure[8]togeth er with our reduced spectrum of HD 52961, on 
which a study (jGielen et al.ll2011ah ha s been published since 
this paper was submitted. Although iGielen et alj (|2011aT ) 
conclude that HR 4049 does not carry the signatures of Ceo 
our analysis indicates that there is evidence of Ceo emission 
as shown in the top-right inset of the figure. Due to some 
contamination of the 17.4 u.m Ceo band from gas phase CO2 
and PAH emission, only an upper limit on the flux in this 
band could be determined. 



6 DISCUSSION 

In this section we discuss the appearance of Ceo in YSOs 
and other objects with particular emphasis on its excita- 
tion (section [671) and formation (section |6.2|I . 



6.1 Excitation of Ceo IR- emission bands 

Two excitation mechanisms for Ceo IR emission with a 
quantitative basis have been discussed in the literature. 
ICami et al.l (|2010l l put forward a model in which thermal 



equilibrium is assumed, which then allows a vibrational tem- 
perature to be deduced. In Tc 1 the temperature was found 
to be unexpectedly low at 330 K and led to the sugges- 
tion that the emitting Ceo molecules are attached to dust 
grains. In t heir study of Ceo em ission from NGC 2023 and 
NGC 7023, ISellgren et ail l|201lD invoked a gas-phase pho- 
toexcitation model commonly applied in interpreting astro- 
nomical infrared emission spectra of PAH molecules. In this 
case absorption of a UV photon is followed by internal re- 
distribution of the energy leading to emission in the mid-IR. 

Both of these quantitative approaches rely on knowl- 
edge of intrinsic infrared band intensities which are 
quite uncertain. F or th e 18.9, 17.4, 8.5 and 7.0 um 
bands ICami et alj l|2010t ) used relative band str e ngths 
of 100 : 48 : 45 : 3 7.8 (iMartin. Koller fc Mihalvl 1 19931 ; 
| Fabianiri996». whereas ISellgren et al. I l|2010fi used those of 
Choi. Kertesz fc Mihalvl feOOOT): IO0T26 : 3 1 : 46. Recently 



Iglesias-Groth. Cataldo fc Manchadd (|201lh measured rela- 
tive band intensities of 100 : 43 : 26.2 : 26.9 at room tem- 
perature in a KBr matrix. The relative band intensities are 
also likely to be temperature-dependent; notably emission 
spectra at thermal equilibrium in the gas phase at ~ 1000 K 
reveal the 18.9 and 17.4 u.m e mission features to have ap- 
proximately the same intensity (|Frum et al.lll99lh . 

Vibrational temperatures for our ob jects have been es - 
timated using the method described by ICami et all (|2010T l 
by plotting ln^u/gu) versus E u /k in an excitation di- 
agram; an example for ISOGAL-P J174639.6-284126 is 
shown in figure [9] Using the flux values given in ta- 
ble [2J which as listed include the [Aril] contribution to the 
7.0 urn band, the Ceo vibrational temperature for ISOGAL- 
P J174639. 6-284126 i s found to b e 670 K (using band 
strengths taken from Cami et alj (|2010|)) or ^ 790 K (us- 
ing ba nd strengths of Iglesias-Groth. Cataldo fe Manchadd 
(|201lf 0. By applying two-component deconvolution fitting 
in 1 RAF, the contribution of [Aril] is estimated to be ~ 
60% and removal leads to lower vibrational temperatures 
of ^ 450 K or ^ 500 K, respectively. Vibrational tempera- 
tures for ISOGAL-P J174639. 6-284126 and the other two 
YSOs are given in table [3] The valu es are in the sam e 
range as determined for PN Tc 1 by ICami et al.l |2010h . 
PNe by iGarcia-Hernandez et ail (|2010T ) and a PPN by 
I Zhang fc Kwokl 1 201 j ). although the low signal-to-noise ra- 
tio of the Ceo bands in SSTGC 372630 renders the appar- 
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Table 3. Observed vi brational band inten sity ratios, inferred vibrational temperatures for YSOs and comparison with predicted photo- 
excited band ratios o flSellgren et alj l|2010h . T, c ih /K and T},~ G /K are the vibrational tem peratures derived from the C60 band strengths 
used bv lCami et al. (2010) and obtained by llglesias-Groth , Cataldo & Manchado (2011), respectively. 



Object! 


^7.04/^18. 9 


^17.4/^18.9 


T c ,/K 

vib' 




ISOGAL-P J174639.6-284126 (CMZ) 


~0.42 a 


0.53 


5g450 a 


^500° 


SSTGC 372630 (CMZ) 


0.70 6 


~0.59 b 


^540 6 


^620 i > 


2MASS J06314796+0419381 


0.29 


0.48 c 


410 


450 


Photon energy/eV 


5 


0.46-0.58 


0.28-0.38 






10 


0.76-0.94 


0.28-0.38 






15 


0.97-1.20 


0.29-0.38 







a Value when 60% contribution to 7.0 |im feature from [Aril] is removed (see text). 

b Silicate and ice absorptions affect continuum level definition. 

c Ratio when contribution of 20% from PAH feature at 17.4 um is removed. 



ently slightly higher temperature values in this YSO less 
certain. We have not conducted our own photexcitation cal- 
culations for infrared emission from Ceo, but predicted in- 
frared band intensity ratios for p hotoexcitation energi es of 
5, 10 and 15 eV are available from lSellgren et"al1 (|201Cf ) and 
are given in table 

We now consider whether either of the thermal or pho- 
toexcitation models hold for these YSOs. Examination of 
table [3] shows that, with one possible exception (/7.04/J18.9 
for SSTGC 372630), the J7.04/J18. 9 and J17.4//18. 9 ratios fall 
outside the ranges calculated by ISellgren et alj l|2010l) for 
photons with 5, 10 and 15 eV. If the C60 molecule were ex- 
cited by UV radiation, the flux ratios derived for ISOGAL-P 
J174639. 6-284126 and 2MASS J06314796+0419381 (see ta- 
ble[3} would suggest the need for photons with e nergies lower 



Excitation Diayi.uii ioi iSOGAL-P J 1746W.6-2K4126 



than 5 eV. A similar result is fou nd for Tc 1 (ICami et al 
2011 ) and for a number of PNe (|Garcfa-Hernandez et al 



20111 ). However, for the PPN IRAS 01005+7910 the band ra- 



tios can be satisfied within this model but with high-e nergy 
photons of around 10 - 15 eV (|Zhang fc Kwokll201ll ). The 
case of the pre-main sequence Herbig Ae/Be star HD 97300 
is notable because even for photon energies as low as 5 
eV, a ratio of 0.46-0.58 for the 7.0/18.9 band ratio is pre- 
dicted. While the 18.9 [im band of C60 is clearly discernible 
in HD 97300, a band at 7.0 [im is not seen which implies 
emission from very cool Ceo- Unfortunately the 17.4 um fea- 
ture is too heavily contaminated by PAH features to allow 
a value for TUb to be deduced. 

From these comparisons, an emission mechanism based 
solely on the absorption of UV radiation followed by re- 
emission in the infrared would seem unlikely for these 
four new Cgo-containing young objects. However, within 
the thermal model, the derived temperatures (table [3} for 
the three YSO objects are not inconsistent with the tem- 
peratu res of up to ^1000 K for warm gas in massive 
YSOs (lLahuis fc van Dishoeckl [ioool ; iBoonman et al.ll2003l ; 
lAn et al.ll201ll ). While the location of the C 6 o in the YSOs 
is not yet established, and there is no evidence for gas-phase 
absorption by HCN, C2H2 or CO2 in the spectra, warm 
gaseous regions could provide a suitable environment for 
thermal excitation of IR emission from Ceo- 

For SSTGC 372630 (CMZ) and 2MASS 
J06314796+0419381 (Rosette nebula) there are S(l) 
and S(2) pure rotational lines of molecular hydrogen. 



1000 1200 1400 1600 1800 

E„/k 



Figure 9. Plot of ln(A u /g u ) versus i? u /k where N u is the 
number of molecules in the upper (u) vibrational level, g u 
is the upper level degeneracy, E u is the l evel energy, and k 
is the Boltzmann constant. For details see ICami et al. I ll2010l) . 
Data for the 18.9, 17.4 and 7.0 ^m bands of C 6 o in the YSO 
ISOGAL-P J174 639.6-284126 are shown. For th i s example the 
band strengths llMartin. Roller fc Mihalvl Il993l ; iFabianl Il996l) 
were used. A vibrational temperature of ~450 K is obtained from 
the slope. 



Assuming that the lines are optically thin and that thermal 
equilibrium holds both for the rotational levels and for 
the ortho/para ratio, then rotational temperatures can be 
deduced. We find T rot values of ~540 and ~370 K for 
SSTGC 372630 and 2MASS J06314796+0419381 which 
are within ~ 100 K of the vibrational temperatures for 
C60 in these objects. In addition to the assumptions listed, 
this difference may have a number of contributing factors 
including the molecular physics of the excitation, the 
allowed vibrational transitions of Ceo and the strongly 
forbidden rotational ones of H2, and uncertainty as to 
whether Ceo and H2 share the same spatial distribution. 
Nevertheless, to our knowledge this is the first time that a 
Ceo vibrational temperature has been compared with the 
internal temperature of another molecule. Observations of 
other (polar) molecules through their rotational spectra in 
these objects would clearly be of interest. 

Further possible Ceo vibrational excitation mecha- 
nisms include shock-induced excitation, possibly during 
Ceo formation from HAC in grain-grain collisions, or 
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through dehydrogenation of HA C while on the grain surface 
i|Garcia-Hernandez et all l201ll ; ICami et all l201ll ). The 
YSOs SSTGC 372630 and 2MASS J06314796+0419381, 
show H2 S(l) and S(2) lines which is an indicator of 
shocked gas, while all three YSOs have strong [Nell], also 
an indicator of shocked gas, but no [Neml emission. In 
their s tudy of C6o-containing PNe, IGarcfa-Hernandez et al.l 
(2011) found very low [Ne in] /[Nell] ratios. This result was 
considered incompatible with the predictions of photoioni- 
sation modelling, but one that could possibly be rationalised 
in terms of shocks. Shocks are expected in YSOs and so 
could provide C60 vibrational excitation. In the case of 
SSTGC 372630, the H 2 lines are shifted by ~200 kms" 1 
relative to the [Nell] line, strongly suggesting a strong 
outflow/shock which is consistent with the likelihood that 
this object is an 'outflow' source as mentioned in section l4~2l 



6.2 Formation of Ceo 

It is well established that C60 is formed i n th e 
later stages of stellar e volut ion (see ICami et al.l l|2010l ). 
IGarcfa-Hernandez et al.l (|2010l ) and reference s in section []]). 
Its p resence in NGC 2023 and NGC 7023 l|Sellgren et al l 
2010) indicates that it might also be formed in situ 
in the ISM, where it ma y also exist in ionised form 
ijFoing fc Ehrenfreundlll994l ). Discovery of C60 in the YSOs 
reported here makes the question as to the mechanism and 
location for C60 formation an even wider one. In this section 
we discuss how C60 could form in situ in pre-main-sequence 
objects or arise from earlier synthesis in post-AGB stars or 
the ISM ( section 16. 2.1[1 . and draw attention to a tentative 
link between C60 and nanodiamonds (section 16. 2. 2[) . 



6.2.1 C 60 m the YSOs and Herbig Ae/Be star HD 97300 

A number of scenarios for the origin of Ceo in YSOs can 
be envisaged. Given that it is clearly formed in post-AGB 
objects, then if the integrity of the Ceo structure were 
maintained on it being expelled to the ISM, it could then 
(re)appear in the spectrum of a YSO following cloud col- 
lapse. For a discussion on the evolution of organic materia l 
through these stages, see lEhrenfreund fc Charnlevl (2000). 
Alternatively Ceo in YS Os could originally have forme d in 
the ISM as discussed by iBettens fc Herbstj jl996l . I1997T ). A 
third scenario is that Ceo forms in situ in YSOs. Recent 
commentaries on the current challenges relating to both 
mechanistic a nd spatial a s pects of Ceo fo rmation have been 
publis hed by ICami et al.l (|201lh and by [G arcf a- Hernandez! 
l|201ll ). Attention is focussed here on in situ formation of 
C 60 in YSOs. 

We consider two main schemes. Firstly, shock- 
induced formation of Ceo might occur in situ in YSOs 
through decomposition of HAC , essentially as described by 
IGarcfa-Hernandez et al.l (2010) for post-AGB objects. How- 
ever, it is not clear whether HAC is present in YSOs. For 
post-AGB objects there is a very strong correlation be- 
tween Ceo and the broad 30 frm feature that is commonly 
attributed to HAC. However, the Ceo-containing reflection 
nebulae NG C 2023 and NGC 7023 do not have the 30 urn 
emission (see lZhang fc Kwokl (|201ll )). and nor do the YSOs 



SSTGC 372630 and ISOG AL-P J174639.6-28 4126, the Her- 
big Ae/Be star HD 97300 jKeller et al.ll200ct) or the mixed- 
chemistry post-AGB objects HR 4049 and HD 52961. Hence, 
if the 30 [im emission does indeed arise from HAC then it is 
hard to see how this particular shock-induced Ceo-formation 
mechanism involving HAC could be operating in young stel- 
lar environment s@ 

Secondly, in most astrophysical sources where Ceo has 
been found, PAH emission is also present. This raises the 
question as to whether there might be a link between PAHs 
(perhaps p articularly in dehydrogenate d form) and fullerene 
formation. iMicelotta. Jones fc Tielensl (|2010 a b) have ar- 
gued from a theoretical standpoint that fullerene formation 
could occur through shock driven processing of PAHs. They 
conclude that shock velocities of 75-100 kms" 1 can signif- 
icantly modify the structure of PAHs, these speeds being 
comparable to those found in massive YSOs. Even at veloc- 
ities as low as 50 kms -1 some carbon atoms are expected 
to be removed from t he carbon framework as could occur 
in C-t y pe shocks (e.g. van den Ancker. Tielens fc Wesseliusl 
l|2000l )). IChuvilin et al.l (|2010l ) have recently shown that un- 



der highly energetic conditions small PAH-like graphene 
sheets can undergo re-arrangement t o form fullerenesLj Al - 
though the laboratory experiments of lChuvilin et all (|2010l ) 
were conducted using transmission electron microscopy in 
which very high energy electrons induce dehydrogenation 
and loss of carbon atoms leading to Ceo formation, high 
energy photons, x-rays and cosmic rays could play an equiv- 
alent role in astrophysical environments. These processes all 
proceed via initial dehydrogenation and loss of carbon atoms 
leading to formation of 5-membered rings. Curvature follows 
with the final stage being closure of the fullerene structure. 
Given that PAHs are present in all of the YSOs, 'top-down' 
mechanisms starting from large PAHs would appear to be 
favoured as formation routes for Ceo in YSOs as well as in 
the Herbig Ae/Be star HD 97300. 



6.2.2 C 60 in post-AGB stars: HR 4049 and HD 52961 

During the search of Spitzer archival data two post-AGB 
sources with Ceo emission at 18.9 u.m were found, HR 4049 
and HD 52961 (see table [1] and figure [8|. These objects 
are unusual in that they have a mixed oxygen-carbon chem- 
istry and so the presence of Ceo is perhaps unexpected. The 
presence of a ban d of Ceo in HP 5296 1 has also been re- 
cently reported by ICielen et al.l (|2011al ) We include these 
objects because of the parallel between the outflows (jets 
and shocks) that are common in the later stages of stel- 
lar evolution and those occurring in young objects. It fol- 
lows that Ceo formation processes in these two very differ- 
ent evolutionary stages may be quite similar. We suggest 



6 The anonymous referee has pointed out that the absence of a 
30 (im feature could be related to the temperature of the bulk of 
the HAC material. 

7 Since this paper w as submitted an arXiv manuscript by 
iBerne fc Ti elens (20l J) has appeared which discusses this work 
with reference to NGC 7023. 

8 Since this paper was submitted an arXiv manuscript by 
iGielen et al.l j2011bT ) has appeared in which a discussion of possi- 
ble formation mechanisms in a mixed oxygen-carbon environment 
is presented. 
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th at 'top-down' formation of Ceo from PAHs as described 
by iMicelotta. Jones fe Tielens! J2010al fbh and seen experi- 
mentallv bv IChuvilin et all (|2Qlof T will likely apply also to 
Ceo formation in mixed chemistry post-AGB objects such 
as HR 4049. 

HR 4049 and HR 52961 exhibit the very rare C- 
H emission features at 3.43/3.53 [ im from H-coate d nan- 
odiamonds llGeballe et al.l 19891: lOudmaiier et al ] 1 19951 : 



iGuillois. Ledoux fc Revnaud 19991 ). Given that they also 
have Ceo emission, this suggests that there may be a connec- 
tion between fullerenes and nanodiamonds. Carbon onions 
are composed of fullerene-like concentric shells and due to 
their similarity with fullerenes may form together in astro- 
physical environments. It has been proposed that astrophys- 
ical carbon onions could act as pressure cells for the for- 
mation of nanodi amonds under highly energetic conditions 
ijGoto et al.ll2009l ). The fact that the only post-AGB objects 
with nanodiamond signatures also have Ceo bands provides 
indirect observa tional support for the nanodiamond forma- 
tion proposal of iGoto et al.l (|2009l ). 



7 SUMMARY 

This paper reports the first detection of the Ceo molecule 
in three young stellar objects and in a Herbig Ae/Be star 
using Spitzer IRS archive observations. Evidence that the 
three objects are YSOs is critically reviewed and their rel- 
ative stage of evolution discussed. Ceo emission in two un- 
usual mixed-chemistry post-AGB objects is also described 
as they share shock-induced processes similar to those active 
in YSOs. A common feature of all of the objects is emission 
from PAHs. 

The origin of the Ceo emission is considered in detail. 
From the measured fluxes in the 7.0, 17.4 and 18.9 (im bands, 
it is found that the band intensities can best be described by 
a thermal rather than a UV photoexcitation model and ap- 
proximate vibrational temperatures were determined. These 
are similar to the rotational temperatures deduced from H2 
line emission found in two of the young objects. 

C60 formation mechanisms are evaluated, focussing par- 
ticularly on photo-induced and shock-induced in situ for- 
mation from PAHs in YSOs. We suggest the most likely 
formation mechanism is through dehydrogenation followed 
by carbon-atom loss from PAHs resulting in five-membered 
ring formation, curvature and closure to make fullerenes. It 
is concluded that shock-induced and UV stimulated decom- 
position are active in generating C60 in the YSOs, with the 
relative contributions varying with evolutionary stage. It is 
thought that shocks play the dominant role in C60 formation 
in the Herbig Ae/Be and post-AGB stars. 
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